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^ ' Abstract We review the experimental and theoretical status of elastic electron scattering and 

I , elastic low-energy photon scattering (with both real and virtual photons) from the nucleon. As 

^ ■ a consequence of new experimental facilities and new theoretical insights, these subjects are 

r—j ' advancing with unprecedented precision. These reactions provide many important insights into 

(~| , the spatial distributions and correlations of quarks in the nucleon. 
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1 General Introduction 

Although nucleons account for nearly all the visible mass in the universe, they 
have a complicated structure that is still incompletely understood. The first indi- 
cation that nucleons have an internal structure, was the 1933 measurement of the 
proton magnetic moment by Frisch & Stern(^. The investigation of the spatial 
structure of the nucleon was initiated by the HEPL (Stanford) experiments in the 
1950s, for which Hofstadter was awarded the 1961 Nobel prize. Several volumes of 
the Annual Review of Nuclear Science^l^ reviewed the status of this field. The 
recent revival of its experimental study through the operational implementation 
of novel instrumentation has instigated a strong theoretical interest. 

Nucleon electro-magnetic form factors (EMFFs) are optimally studied through 
the exchange of a virtual photon, in elastic electron-nucleon scattering. The 
momentum transfer to the nucleon can be selected to probe different scales of the 
nucleon, from integral properties such as the charge radius to scaling properties 
of its internal constituents. Polarization instrumentation, polarized beams and 
targets, and the measurement of the recoil polarization have been essential in the 
accurate separation of the charge and magnetic form factors and in studies of the 
elusive neutron charge form factor. 

Exclusive Compton scattering on a nucleon refers to the reactions l'^', 
where either photon may be real or virtual. In general, the Compton amplitude 
depends on the full complexity of the dynamics of the excitation spectrum of 
the nucleon. However, in a number of special kinematic domains, observables 
with a particularly simple interpretation can be extracted from the Compton 
amplitude. In this review, we present the experimental and theoretical status of 
real (RCS) and virtual (VCS) Compton scattering for the study of generalized 
polarizabilities, which measure the spatial distribution of the response of the 
nucleon to external electromagnetic fields. A thorough discussion of the rapid 
developments in high energy Compton scattering, in both the deep virtual and 
hard scattering limits, is beyond the scope of this review. 

2 Nucleon Form Factors 

2.1 Theory of Electron Scattering and Form Factor Measure- 



The nucleon EMFFs are of fundamental importance for the understanding of the 
nucleon's internal structure. Under Lorentz invariance, spatial symmetries, and 
charge conservation, the most general form of the electromagnetic current inside 
a nucleon can be written as: 



where -Fi denotes the helicity non-flip Dirac form factor, F2 the helicity flip Pauli 
form factor, = —q^, and k the nucleon anomalous magnetic moment. The 
remaining variables are defined in Figure ^ The second term, usually referred 
to as the Foldy contribution, is of relativistic origin. It is useful to introduce the 
isospin form-factor components, corresponding to the isoscalar (s) and isovector 
{v) response of the nucleon. 
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Fl = ^{F!' + Fn; F^ = \{Ff-F^), (z = l,2). (2) 




Figure 1: The Feynman diagram for the scattering of an electron with four- 
momentum p = (Ee,p) through an angle 9e off a nucleon with mass Mjy and 
four-momentum P. In this diagram a single virtual photon with four-momentum 
q = p — p' = (lo, q) is exchanged. The four-momenta of the scattered electron and 
nucleon are p' = {E'^,p') and P' , respectively. 

The form factors can be continued analytically into the complex plane and can 
be related in different regions through a dispersion relation of the form 

^w^irfe.., (3) 

vr J to V —t 

with t = —Q^, to = 9(4)M^ for the isoscalar (isovector) case and the pion 
mass. In the positive Q^-region, called spacelike, form factors can be measured 
through electron scattering, in the negative Q'^-region, called timelike, form fac- 
tors can only be measured through the creation or annihilation of a A^iV-pair. 
In plane-wave Born approximation, the cross section for elastic electron-nucleon 
scattering can be expressed in the Rosenbluth@ formula as: 

^ = aM[{Fl + + 2r(Fi + nF.f tan2(|)], (4) 

where r = QV(4M^) and cjm = is the Mott cross section for 

scattering off a point-like particle, with uqed denoting the fine-structure con- 
stant. The remaining variables are defined in Figure ^ IIofstadter(^ determined 
the values of Fi and F2 from measurements at different scattering angles but at 
the same values of by drawing intersecting ellipses. Hand, Miller and Wil- 
son® pointed out that a simple algebraic separation is possible if one expresses 
the Rosenbuth formula in an alternate form: 



(5) 

with e = 1/[1 -|- 2(1 -|- T)tan^(y)] the linear polarization of the virtual photon. 
They further noted that Ge and Gm were identical to the electric and magnetic 
form factors, discussed earlier by Ernst, Sachs and Wali(|7|): 
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GM{Q^)=Fi{Q^) + nF2{Q^); G^,f (0) = (6) 

with ^p^n denoting the magnetic moment of the proton and neutron, respectively. 
Equation [S] ihustrates that and G^^ can be determined separately by per- 
forming cross-section measurements at fixed as a function of e, over a range of 
{9e,Ee) combinations (Rosenbluth separation). Hand, Miller and Wilson further 
noted that in the Breit frame, which for elastic scattering is equivalent to the 
centre-of-mass frame, the electromagnetic current of the proton simplifies into 
the expression: 

J^Af = GEiQ^h^ + GM{Q^)icT^''qu. (7) 

In this reference frame the Sachs form factors can be identified with the Fourier 
transform of the nucleon charge and magnetization density distributions. 

Through the mid-1990s practically all available proton EMFF data had been 
collected using the Rosenbluth separation technique. This experimental proce- 
dure requires an accurate knowledge of the electron energy and the total lumi- 
nosity. In addition, because the contribution to the elastic cross section is 
weighted with Q^, data on G^ suffer from increasing systematic uncertainties 
with increasing Q^-values. The then available world data setl'^ was compared 
to the so-called dipole parametrization Gd, which corresponds to two poles with 
opposite sign close to each other in the time-like region. In coordinate space Gd 
corresponds to exponentially decreasing radial charge and magnetization densi- 
ties, albeit with a non-physical discontinuity at the origin: 

/ A2 \' 

Gd= { To 7^ with A = 0.84 GeV and Q in GeV. (8) 

\^ + (y^ y 

For G^, G\j/fi.p and G\[/fin the available data agreed to within 20% with 
the dipole parametrization. Both the G^; and the G^/fj,p data could be fit- 
ted adequately with an identical parametrization. However, the limitation of 
the Rosenbluth separation was evident from the fact that different data sets for 
^pG^/G^ scattered by up to 50% at Q^-values larger than 1 GeV^ (Figure IS}. 
Although no fundamental reason has been found for the success of the dipole 
parametrization, it is still used as a base line for comparison of data because it 
takes out the largest variation with and enables smaller differences to be seen. 



2.2 Instrumentation for Form Factor Measurements 

More than 40 years ago Akhiezer et al.pSj) (followed 20 years later by Arnold 
et al.(16 )) showed that the accuracy of nucleon charge form-factor measurements 
could be increased significantly by scattering polarized electrons off a polarized 
target (or equivalently by measuring the polarization of the recoiling proton). 
However, it took several decades before technology had sufficiently advanced to 
make the first of such measurements feasible, and only in the past few years 
has a large number of new data with a significantly improved accuracy become 
available. The next few sections introduce the various techniques. For G^ mea- 
surements, the highest figure of merit at Q^-values larger than a few GeV^ is 
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Figure 2: The ratio iipG^^/G\^ from Rosenbluth separation. Data are from Ref- 
erences (|ni llUl[TTlll21ll31ll4|) . The errors shown are the quadratic sum of the 
statistical and systematic contributions. 



obtained with a focal plane polarimeter. Here, the Jacobian focusing of the re- 
coiling proton kinematics allows one to couple a standard magnetic spectrometer 
for the proton detection to a large-acceptance non-magnetic detector for the de- 
tection of the scattered electron. For studies of one needs to use a magnetic 
spectrometer to detect the scattered electron in order to cleanly identify the re- 
action channel. As a consequence, the figure of merit of a polarized ^He target is 
comparable to that of a neutron polarimeter. 

2.2.1 Polarized Beam Various techniques are available to produce polar- 
ized electron beams, but photo-emission from GaAs has until now proven to be 
optimal(^7J. A thin layer of GaAs is illuminated by a circularly polarized laser 
beam of high intensity, which preferentially excites electrons of one helicity state 
to the conductance band through optical pumping. The helicity sign of the laser 
beam can be flipped at a rate of tens of Hertz by changing the high voltage on a 
Pockels cell. The polarized electrons that diffuse to the photocathode surface are 
then extracted by a 50-100 kV potential. An ultra-high vacuum environment is 
required to minimize surface degradation of the GaAs crystal by backstreaming 
ions. Initially, the use of bulk GaAs limited the maximum polarization to 50% 
because of the degeneracy of the P3/2 sublevels. This degeneracy is removed by 
introducing a strain in a thin layer of GaAs deposited onto a thicker layer with 
a slightly different lattice spacing. Although such strained GaAs cathodes have 
a significantly lower quantum efficiency than bulk GaAs cathodes, this has been 
compensated by the development of high- intensity diode or Ti-sapphire lasers. 
Polarized electron beams are now reliably available with a polarization close to 
80% at currents of > 100 //A. 
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The polarized electrons extracted from the GaAs surface are first pre- accelerated 
and longitudinally bunched and then injected into an accelerator. Typically, the 
polarization vector of the electrons is manipulated in a Wien filter so that the 
electrons are fully longitudinally polarized at the target. If the beam is injected 
into a storage ring for use with an internal target, a Siberian snake(|18)l is needed 
to compensate for the precession of the polarization. 

Three processes are used to measure the beam polarization: MottfTO') scat- 
tering, Moller (l2n|) scattering or Compton(21| scattering. Any of these results 
in a polarimeter with an accuracy approaching 1%. In a Mott polarimeter the 
beam helicity asymmetry is measured in scattering polarized electrons off atomic 
nuclei. This technique is limited to electron energies below ~ 20 MeV and mul- 
tiple scattering effects have to be estimated by taking measurements at different 
target foil thicknesses. In a M0ller polarimeter polarized electrons are scattered 
off polarized atomic electrons in a magnetized iron foil. In this technique the 
major uncertainties are in the corrections for atomic screening and in the foil 
magnetization, unless the polarizing field is strong enough to saturate the mag- 
netization. A potentially superior alternative()22|) has been proposed in which the 
electrons are scattered off a sample of atomic hydrogen, polarized to a very high 
degree in an atomic beam, and trapped in a superconducting solenoid. Finally, 
in a Compton polarimeter the beam helicity asymmetry is measured in scattering 
polarized electrons off an intense beam of circularly polarized light, produced by 
trapping a laser beam in a high-finesse Fabry-Perot cavity. The electron beam in 
a storage ring is sufficiently intense that a laser beam can be directly scattered 
off the electron beam without the use of an amplifying cavity. Only the last 
two methods, the atomic hydrogen M0ller and the Compton polarimeter, have 
no effect on the quality of the electron beam and thus can be used continuously 
during an experiment. 

2.2.2 Polarized Targets In polarized targets for protons two different 
techniques are used, depending on the intensity of the electron beam. In storage 
rings where the circulating beam can have an intensity of 100 mA or more, but the 
material interfering with the beam has to be minimized, gaseous targets are used, 
whereas in external targets solid targets can be used. Because free neutrons are 
not available in sufficient quantity, effective targets, such as deuterium or ^He, 
are necessary, and the techniques used to polarize the deuteron are similar to 
those used for the proton. For ^He gaseous targets are used both in internal and 
external targets. 

Solid polarized targets that can withstand electron beams with an intensity of 
up to 100 nA all use the dynamic nuclear polarization technique(23 ). A hydroge- 
nous compound, such as NH3 or LiD, is doped, e.g. by radiation damage, with 
a small concentration of free radicals. Because the occupation of the magnetic 
substates in the radicals follows the Boltzmann distribution, the free electrons 
are polarized to more than 99% in a ~ 5 T magnetic field and at a ~ 1 K tem- 
perature. A radiofrequency (RF) field is then applied to induce transitions to 
states with a preferred orientation of the nuclear spin. Because the relaxation 
time of the electrons is much shorter than that of the nuclei, polarized nuclei 
are accumulated. This technique has been successful in numerous deep-inelastic 
lepton scattering and nucleon form-factor experiments; it has provided polarized 
hydrogen or deuterium targets with an average polarization of ~ 80% or ~ 30%, 
respectively. 
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Internal hydrogen/deuterium targets (|24() are polarized by the atomic beam 
source (ABS) technique, which relies on Stern-Gerlach separation and RF transi- 
tions. First, a beam of atoms is produced in an RF dissociator through a nozzle 
cooled with liquid nitrogen. Then, atoms with different electron spin direction are 
separated through a series of permanent (or superconducting) sextupole magnets 
and transitions between different hyperfine states are induced by a variety of RF 
units. The result is a highly polarized beam with a flux up to 10^"^ atoms/s. This 
beam is then fed into an open-ended storage cell, which is cooled and coated to 
minimize recombination of the atoms bouncing off the cell walls. The circulating 
electron beam, passing through the long axis of the storage cell, encounters only 
the flowing atoms. The polarization vector is oriented with a set of coils, produc- 
ing a field of ~ 0.3 T in order to minimize depolarization by the RF structure 
of the circulating electron beam. The diameter of the storage cell is determined 
by the halo of the electron beam. A target thickness of 2 x 10^^ nuclei/cm^ has 
been obtained at a vector polarization of more than 80%. 

Polarized hydrogen or deuterium atoms can also be produced by spin-exchange 
collisions between such atoms and a small admixture of alkali atoms that have 
been polarized by optical pumping. The nucleus is then polarized in spin- 
temperature equilibrium. Although the nuclear polarization obtained in such 
a laser driven source (LDS) is smaller than through the ABS technique, the flux 
can be more than 10^® atoms/s. Moreover, an LDS offers a more compact design 
than an ABS. A figure of merit comparable to that of the ABS at the HERMES 
experiment has recently been achieved by the MIT group( 25). 

Polarized ^He is attractive as an effective polarized neutron target because 
its ground state is dominated by a spatially symmetric s-state in which the 
proton spins cancel, so that the spin of the ^He nucleus is mainly determined 
by that of the neutron. Corrections for the (small) d-state component and for 
charge-exchange contributions from the protons can be calculated accurately at 
Q^-values smaller than 0.5 GeV'^^^ and larger than ~ 2 GeV^(j22j). Direct 
optical pumping of ^He atoms is not possible because of the energy difference 
between the ground state and the first excited state. Instead ^He is polar- 
ized, either by first exciting the atoms to a metastable 2^Si state and optically 
pumping that state, which then transfers its polarization to the ground state by 
metastability-exchange collisions, or by optically pumping a small admixture of 
rubidium atoms, which then transfer their polarization to the ^He atoms through 
spin-exchange collisions. In internal targets only the metastability-exchange tech- 
nique has been used because of the possible detrimental effects of the rubidium 
admixture on the storage ring environment. With beam on target, polarization 
values of up to 46% at target thicknesses of 1 x 10^^ nuclei/cm^ have been ob- 
tained. For external targets the spin-exchange technique (j2H|) has been used to 
optically pump a glass target cell filled with 10 atm of ^He with a 0.1% rubidium 
admixture. After the spin-exchange collisions the polarized ^He diffuses into a 25 
cm long cell which the electron beam traverses. Polarizations in excess of 40% 
have been reached with beam on target. A pair of 5 mT Helmholz coils is used to 
orient the polarization vector, and care must be taken to minimize depolarizing 
magnetic field gradients. Alternatively, the metastability techniaue H29() has been 
used to polarize ^He under atmospheric pressure which is then compressed to a 
density of more than 6 atm. 
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2.2.3 Recoil Polarimeters Focal-plane polarimeters have long been used 
at proton scattering facilities to measure the polarization of the scattered pro- 
ton. In such an instrument ()28() the azimuthal angular distribution is measured of 
protons scattered in the focal plane of a magnetic spectrometer by an analyzer, 
which often consists of carbon. From this angular distribution the two polar- 
ization components transverse to the proton momentum can be derived. The 
analyzer is preceded by two detectors, most often wire or straw chambers, to 
measure the track of the incident proton; it is followed by two more detectors to 
track the scattered particle. The thickness of the analyzer is adjusted to the pro- 
ton momentum, limiting multiple scattering while optimizing the figure of merit. 
In order to determine the two polarization components in the scattering plane at 
the target, care must be taken to accurately calculate on an event-by-event basis 
the precession of the proton spin in the magnetic field of the spectrometer. 

Neutron polarimeters follow the same basic principle. Here, plastic scintilla- 
tor material is used as an active analyzer, preceded by a veto counter to discard 
charged particles. This eliminates the need for the front detectors. Sets of scin- 
tillator detectors are used to measure an up-down asymmetry in the scattered 
neutrons, which is sensitive to a polarization component in the scattering plane, 
perpendicular to the neutron momentum. In modern neutron polarimeters (jSUj) 
the analyzer is preceded by a dipole magnet, with which the neutron spin can be 
precessed. 

2.3 Experimental Results 

2.3.1 Proton Electric Form Factor In elastic electron-proton scatter- 
ing a longitudinally polarized electron will transfer its polarization to the recoil 
proton. In the one-photon exchange approximation the proton can attain only 
polarization components in the scattering plane, parallel (Pi) and transverse (Pj) 
to its momentum. This can be immediately seen from the expression for the 
proton current in the Breit frame, which separates into components proportional 
to G^; and Gm to which reference was made in eq. The ratio of the charge and 
magnetic form factors is directly proportional to the ratio of these polarization 
components(|,31 J : 



The polarization-transfer technique was used for the first time by Milbrath 
et al. 1)32(1 at the MIT-Bates facility. The proton form factor ratio was measured at 
Q^-values of 0.38 and 0.50 GeV^ by scattering a 580 MeV electron beam polarized 
to ~ 30%. A follow-up measurement was performed at the MAMI facility (jHSJ) at 
a Q2.value of 0.4 GeV^. 

The greatest impact of the polarization-transfer technique was made by the 
two recent experimentsf34''35') in Hall A at Jefferson Lab, which measured the ra- 
tio G^/G^j in a Q^-range from 0.5 to 5.6 GeV^. Elastic ep events were selected by 
detecting electrons and protons in coincidence in the two identical high-resolution 
spectrometers. At the four highest Q^-values a lead-glass calorimeter was used to 
detect the scattered electrons in order to match the proton angular acceptance. 
The polarization of the recoiling proton was determined with a focal-plane po- 
larimeter in the hadron spectrometer, consisting of two pairs of straw chambers 
with a carbon or polyethylene analyzer in between. The data were analyzed in 




Pt Ee + E'^ 

Pi 2M 




(9) 
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Figure 3: The ratio /UpG^/G^j from polarization transfer . com- 
pared to recent Rosenbluth data (|39[l37j) and the reanalysis by Arrington(38) of 
older SLAC data. 



bins of each of the target coordinates. No dependence on any of these variables 
was observed(34). Figure shows the results for the ratio UpG^/G^ . The most 
striking feature of the data is the sharp, practically linear decline as increases: 

Mp^^ = 1 - 0.13(Q2 _ 0.29GeV2). (10) 

Since it is known that G^f closely follows the dipole parametrization, it fol- 
lows that G^ falls more rapidly with than Go- This significant fall-off of the 
form-factor ratio is in clear disagreement with the results from the Rosenbluth 
extraction. Arrington(|38|) has performed a careful reanalysis of earlier Rosen- 
bluth data. He selected only experiments in which an adequate e-range was 
covered with the same detector. The results (Figure do not show the large 
scatter seen in Figure l2j Recently, Christy et al. (|37|) analyzed an extensive data 
set on elastic electron-proton scattering collected in Hall C at Jefferson Lab as 
part of experiment E99-119. The results are evidently in good agreement with 
Arrington's reanalysis. Qattan et al. (39) performed a high-precision Rosenbluth 
extraction in Hall A at Jefferson Lab, designed specifically to significantly reduce 
the systematic errors compared to earlier Rosenbluth measurements. The main 
improvement came from detecting the recoiling protons instead of the scattered 
electrons, so that the proton momentum and the cross section remain practically 
constant when one varies e at a constant Q^-value. In addition, possible depen- 
dences on the beam current are minimized. Special care was taken in surveying 
the angular setting of the identical spectrometer pair. One of the spectrometers 
was used as a luminosity monitor during an e scan. The results(39) of this ex- 
periment, covering Q^-values from 2.6 to 4.1 GeV^, are in excellent agreement 
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with previous Rosenbluth results. This basically rules out the possibility that 
the disagreement between Rosenbluth and polarization-transfer measurements of 
the ratio G^/G^^ is due to an underestimate of e-dependent uncertainties in the 
Rosenbluth measurements. 

2.3.2 Two-Photon Exchange In order to resolve the discrepancy be- 
tween the results for G^/G^ from the two experimental techniques, an e-dependent 
modification of the cross section is necessary. In two-(or more-)photon exchanges 
(TPE) the nucleon undergoes a first virtual photon exchange which can lead to 
an intermediate excited state and then a second one or more, finally ending back 
in its ground state (Figure 2)) • The TPE contributions to elastic electron scat- 
tering have been investigated both experimentally and theoretically for the past 
fifty years. In the early days such contributions were called dispersive effects IIQj. 
Lately, they have been relocated to radiative corrections in the so-called box dia- 
gram. Almost all analyses with the Rosenbluth technique have applied radiative 
corrections using the formulae derived by Mo & Tsai (|41() that only include the 
infrared divergent parts of the box diagram (in which one of the two exchanged 
photons is soft). Thus, terms in which both photons are hard (and which depend 
on the hadronic structure) have been ignored. 



P P' 




P' 



Figure 4: The Feynman diagrams depicting two-photon exchanges. 

The most stringent tests of TPE on the nucleon have been carried out by mea- 
suring the ratio of electron and positron elastic scattering off a proton. Correc- 
tions due to TPE will have a different sign in these two reactions. Unfortunately, 
this (e'''e~) data set is quite limited (|44)1 . only extending (with poor statistics) 
up to a Q^-value of ~ 5 GeV^, whereas at Q^-values larger than ~ 2 GeV^ ba- 
sically all data have been measured at e-values larger than ~ 0.85. Other tests, 
also inconclusive, searched for non-linearities in the e-dependence or measured 
the transverse (out-of-plane) polarization component of the recoiling proton, of 
which a non-zero value would be a direct measure of the imaginary part of the 
TPE amplitude. 

Several studies have provided estimates of the size of the e-dependent cor- 
rections necessary to resolve the discrepancy. Because the fall-off of the form- 
factor ratio is linear with Q^, and the Rosenbluth formula also shows a linear 
dependence of the form-factor ratio (squared) with through the r-term, a Q^- 
independent correction linear in e would cancel the disagreement. An additional 
constraint that any e-dependent modification must satisfy, is the (e"'"e~) data set. 
Guichon & Vander haeghen H42() introduced a general form of a TPE contribution 
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from the so-called box diagram in radiative corrections into the amplitude for 
elastic electron-proton scattering. This resulted in the following modification of 
the Rosenbluth expression: 

~ 2 

daocT + e^ + 2e{T + ^)Y2^, (11) 
Gm 

where Y2'y = , ^fj^- and Gm, F2 and F3 are equal to , F2 and 0, respectively, 
in the Born approximation. and the "two-photon" form factors Ge and 

Gm were fitted (|42|1 to the Rosenbluth and polarization transfer data sets. This 
resulted in a value of ^ 0.03 for ^27 with very little e- or Q^-dependence. 
ArringtonHSD performed a fit to the complete data set, investigating two different 
modifications to the cross section with a Q^-independent linear e-dependence of 6 
% over the full e-range. Both modifications have the same e-dependence, but one 
does not modify the cross section at small values of e, whereas the other leaves 
the cross section unchanged at large values of e. He found that the second gave 
a much better description of the complete data set. Moreover, it was in good 
agreement with the data set for the ratio of electron-proton and positron-proton 
elastic scattering. 

Blunden et al. (j45(l carried out the first calculation of the elastic contribution from 
TPE effects, albeit with a simple monopole Q^-dependence of the hadronic form 
factors: G{Q'^) = A^/(Q^ + A^). They obtained a practically Q^-independent 
correction factor with a linear e-dependence that vanishes at forward angles (e = 
1). However, the size of the correction only resolves about half of the discrepancy. 
A later calculation(W. Melnitchouk, private communication) which used a more 
realistic form factor behavior, resolved up to 80% of the discrepancy. 
A different approach was used by Chen et al. f?6|) , who related the elastic electron- 
nucleon scattering to the scattering off a parton in a nucleon through generalized 
parton distributions. TPE effects in the lepton-quark scattering process are cal- 
culated in the hard-scattering amplitudes. The handbag formalism of the gener- 
alized parton distributions is extended in an unfactorized framework in which the 
x-dependence is retained in the scattering amplitude. Finally, a valence model is 
used for the generalized parton distributions. The results for the TPE contribu- 
tion fully reconcile the Rosenbluth and the polarization-transfer data and retain 
agreement with positron-scattering data. 

Hence, it is becoming more and more likely that TPE processes have to be 
taken into account in the analysis of Rosenbluth data and that they will af- 
fect polarization-transfer data only at the few percent level. Of course, further 
effort is needed to investigate the model-dependence of the TPE calculations. Ex- 
perimental confirmation of TPE effects will be difficult, but certainly should be 
continued. The most direct test would be a measurement of the positron-proton 
and electron-proton scattering cross-section ratio at small e- values and Q^-values 
above 2 GeV^. Positron beams available at storage rings are too low in either 
energy or intensity, but a measurement in the CLAS detector at Jefferson Lab, 
a more promising venue, has been proposed|42lj). A measurement of the beam 
or target single-spin asymmetry normal to the scattering plane, which directly 
accesses the imaginary part of the box diagrams, would provide a sensitive test 
of TPE calculations. Also, real and virtual Compton scattering data can provide 
additional constraints on calculations of TPE effects in elastic scattering. Rosen- 
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bluth analyses have so far been restricted to simple PWBA, Coulomb distortion 
effects should certainly be included too. Additional efforts should be extended 
to studies of TPE effects in other longitudinal-transverse separations, such as 
proton knock-out and deep-inelastic scattering (DIS) experiments. 
2.3.3 Proton Magnetic Form Factor An extensive data set(I3HI) with a 
good accuracy is available up to a Q^-value of more than 30 GeV^ from unpolar- 
ized cross-section measurements (Figure IS)). Because dominates in a Rosen- 
bluth extraction at larger Q^-values, the G^j data have only a minor sensitivity to 
the discrepancy between the Rosenbluth extraction and the polarization-transfer 
technique. Brash et al.lEI have shown that the G^^ data must be renormalized 
upwards by ~ 2% if one assumes the polarization-transfer data to be correct. 
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Figure 5: The proton magnetic form factor G^^, in units of fJ-pGo, as a function 
of Q"^. Data are from References (|i48..49..50..13..14j) . 



2.3.4 Neutron Magnetic Form Factor Early data on G%j were ex- 
tracted from inclusive quasi-elastic scattering off the deuteron. However, model- 
ing of the deuteron wave function, required to subtract the contribution from the 
proton, resulted in sizable systematic uncertainties. A significant break-through 
was made by measuring the ratio of quasi-elastic neutron and proton knock-out 
from a deuterium target. This method has little sensitivity to nuclear binding 
effects and to fluctuations in the luminosity and detector acceptance. The basic 
set-up used in all such measurements is very similar: the electron is detected in a 
magnetic spectrometer with coincident neutron/proton detection in a large scin- 
tillator array. The main technical difficulty in such a ratio measurement is the 
absolute determination of the neutron detection efficiency. Such measurements 
have been pioneered for Q^-values smaller than 1 GeV^ at Mainz ()531I541I55|) and 
Bonn H56|l . The Mainz data are 8%-10% lower than those from Bonn, at 
variance with the quoted uncertainty of ~2%. This discrepancy would require a 



Nucleon Form Factors and Conipton Scattering 



13 



16%-20% error in the detector efficiency. 

A study of at Q^-values up to 5 GeV^ has recently been completed in Hall B 
by measuring the neutron/proton quasi-elastic cross-section ratio using the CLAS 
detector (157)) ■ A hydrogen target was in the beam simultaneously with the deu- 
terium target. This made it possible to measure the neutron detection efficiency 
by tagging neutrons in exclusive reactions on the hydrogen target. Preliminary 
results (|57)l indicate that is within 10% of Gd over the full Q^-range of the 
experiment (0.5-4.8 GeV^). 

Inclusive quasi-elastic scattering of polarized electrons off a polarized ^He target 
offers an alternative method to determine G^.j through a measurement of the 
beam asvmmetrv(|58|) 



where 6* and (j)* are the polar and azimuthal target spin angles with respect 
to q, Ri denote various nucleon response functions, and Vi the corresponding 
kinematic factors. By orienting the target polarization parallel to q, one measure 
Rt/, which in quasi-elastic kinematics is dominantly sensitive to (GJ^)^. For the 
extraction of Gjy^ corrections for the nuclear medium i2.6| are necessary to take 
into account effects of final-state interactions and meson-exchange currents. The 
first such measurement was carried out at Bates|59). Recently, this technique 
was used to measure G^j in Hall A at Jefferson Lab in a Q'^-range from 0.1 to 
0.6 GeV^ |60{) . This experiment provided an independent, accurate measurement 
of Gjyj at Q^-values of 0.1 and 0.2 GeV^, in excellent agreement with the Mainz 
data. At the higher Q^-values G^j could be extracted (j6ip in plane wave impulse 
approximation, since final-state interaction effects are expected to decrease with 
increasing Q^. 

Figure shows the results of all completed G^ experiments. 
2.3.5 Neutron Electric Form Factor Analogously to G^j, early G^- 
experiments used (quasi-)elastic scattering off the deuteron to extract the lon- 
gitudinal deuteron response function. Due to the smallness of G^, the use of 
different nucleon-nucleon potentials resulted in a 100% spread in the resulting 
G^ values (|65]). In the past decade a series of double-polarization measurements 
of neutron knock-out from a polarized ^H or '^He target have provided accurate 
data on G^ . The ratio of the beam-target asymmetry with the target polarization 
perpendicular and parallel to the momentum transfer is directly proportional to 
the ratio of the electric and magnetic form factors, 

where Px and denote the polarization component perpendicular and parallel 
to q. A similar result is obtained with an unpolarized deuteron target when one 
measures the polarization of the knocked-out neutron as a function of the angle 
over which the neutron spin is precessed with a dipole magnet: 



here, 6 denotes the precession angle where the measured asymmetry is zero. 



A = - 



(cos 0*vt'Rt' + 2 sin 6** cos 9*vtl'Rtl') 
vlRl + vtRt 



(12) 




(14) 
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Figure 6: The neutron magnetic form factor G^j, in units of UnGo, as a function 
of Q^. Results from ^He are indicated by open symbols. Data are from References 
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Again, the first such measurements were carried out at Bates, both with a polar- 
ized ^He target and with a neutron polarimeter. Figure [3 shows results obtained 

through all three reactions ■^H(e, e'n), ^H(e, e'n) and ^He(e, e'n). At low Q'^-values 
corrections for nuclear medium and rescattering effects can be sizeable: 65% for 
at 0.15 GeV^ and 50% for ^He at 0.35 GeV^. These corrections are expected 
to decrease significantly with increasing Q, although no reliable calculations are 
presently available for ^He above 0.5 GeV^. There is excellent agreement between 
the results from the different techniques. Moreover, medium effects have clearly 
become negligible at ~ 0.7 GeV'^, even for ^He. The latest data from Hall C at 
Jefferson Lab, using either a polarimeter or a polarized target (|69 ( I71 |1 . extend 
up to ~ 1.5 GeV^ with an overall accuracy of ~10%, in mutual agreement. 
From ~ 1 GeV^ onwards appears to exhibit a Q^-behavior similar to that of 
G^. Schiavilla &: Sick ((75|l have extracted from available data on the deuteron 
quadrupole form factor i^c2(Q^) with a much smaller sensitivity to the nucleon- 
nucleon potential than from inclusive (quasi-)elastic scattering. The 30-years-old 
Galster Darametrization ((76,) continues to provide a fortuitously good description 
of the data. 

2.3.6 Timelike Form Factors In the timelike region EMFF measure- 
ments have been made at electron-positron storage rings or by studying the in- 
verse reaction (only for the proton form factors), antiproton annihilation on a 
hydrogen target. The rather limited data set on timelike form factors is shown in 
Figure |B1 The quality of the data does not allow a separation of the charge and 
magnetic form factors; Gm has been extracted from the data using the G£;-values 
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Figure 7: The neutron electric form factor as a function of Q^. Results from 

^He are indicated by open symbols. Data are from References H66ll67llF8ll3UllH^ 
|701[7Il|711I731lini|711). The full curve shows the Galster dTHJ) parametrization; the 
dashed curve represents the Q^-behavior of G^. 



calculated by lachello & WanfTT'). Clearly Gd which gives a very good descrip- 
tion of the spacelike magnetic form factors, does not describe the data in the 
timelike region, at least from threshold down to -6 GeV^. lachello & Wan H77|) . 
Hammer et al. (|78() and Dubnicka et al.(79) have carried out an analytic con- 
tinuation of their VMD calculations (section 12. 4|) . lachello's model provides a 
consistent description of the magnetic form factors in the timelike region. An ex- 
tension of the data set in the timelike region and of theoretical efforts to obtain 
a consistent description of all EMFFs in both the space- and timelike regions is 
highly desirable. 

2.3.7 Experimental Review and Outlook In recent years highly ac- 
curate data on the nucleon EMFFs have become available from various facilities 
around the world, made possible by the development of high luminosity and novel 
polarization techniques. These have established some general trends in the Q^- 
behavior of the four EMFFs. The two magnetic form factors G^,j and G^^ are 
close to identical, following G^i to within 10% at least up to 5 GeV^, with a 
shallow minimum at ~ 0.25 GeV^ and crossing Gd at ~ 0.7 GeV^. G^/G^j 
drops linearly with and G^ appears to drop from ~ 1 GeV^ onwards at 
the same rate as G^. Measurements that extend to higher Q^-values and offer 
improved accuracy at lower Q^-values, will become available in the near future. 
In Hall C at Jefferson Lab Perdrisat et al.fSO) will extend the measurements 
of G^/G^ to 9 GeV^ with a new polarimeter and large-acceptance lead-glass 
calorimeter. Wojtsekhowski et al. (|81|) will measure in Hall A at Q^-values 
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Figure 8: The magnetic form factors (divided by Go) in the time-like region as 
a function of Q^, compared to the calculations by Iachello (|77() . Hammer ((751) 
Dubnicka (|79|) ■ See Reference f(77j) for the references to the experimental data. 



of 2.4 and 3.4 GeV^ using the ^He(e, e'n) reaction with a 100 msr electron spec- 
trometer. The Bates Large Acceptance Spectrometer Toroid facility (BLAST, 
http://www.mitbates.mit.edu) at MIT with a polarized hydrogen and deuteron 
target internal to a storage ring will provide highly accurate data on and 
in a Q^-range from 0.1 to 0.8 GeV^. Gao et al.(jES| have shown that the proton 
charge radius can be measured with unprecedented precision by measuring the 
ratio of asymmetries in the two sectors of the BLAST detector. Thus, within a 
couple of years data with an accuracy of 10% or better will be available up to 
a Q^-value of 3.4 GeV^. Once the upgrade to 12 GeV H83() has been implemented 
at Jefferson Lab, it will be possible to extend the data set on and G^ to 14 
GeV2 and on G| to 8 GeV^. 

The charge and magnetization rms radii are related to the slope of the form factor 
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< r|; >= J p{r)r^dr 

<r|,>=/Mr)r** = 

with p(r) (^(r)) denoting the radial charge (magnetization) distribution. Table 
n lists the results. For an accurate extraction of the radius Sick(84) has shown 
that it is necessary to take into account Coulomb distortion effects and higher 
moments of the radial distribution. His result for the proton charge radius is 
in excellent agreement with the most recent three- loop QED calculation (jSSj) of 
the hydrogen Lamb shift. Within error bars the rms radii for the proton charge 
and magnetization distribution and for the neutron magnetization distribution 
are equal. The value for the neutron charge radius was obtained by measuring 
the transmission of low-energy neutrons through liquid 208p^ g,^^^ ^^^^Bi. The 
Foldy term = —0.126 fm^ is close to the value of the neutron charge radius. 
Isgur(87) showed that the Foldy term is canceled by a first-order relativistic 
correction, which implies that the measured value of the neutron charge radius 
is indeed dominated by its internal structure. 



-6 



6 dG(Q2) 



fi dQ^ 



Q2=0 



Q2=o 



(15) 



Table 1: Values for the nucleon charge and magnetization radii 



Observable value ib error Reference 





>l/2 


0.895 ± 0.018 fm 


(|g4|) 


< K,? 


>l/2 


0.855 ± 0.035 fm 


(jHl 


<{rl) 


2> 


- 0.119 ± 0.003 fm^ 


(uni) 


< {rli? 


>l/2 


0.87 ± 0.01 fm 


(55) 



In the Breit frame the nucleon form factors can be written as Fourier transforms 
of their charge and magnetization distributions. However, if the wavelength of 
the probe is larger than the Compton wavelength of the nucleon, i.e. if \Q\ > Mj\f, 
the form factors are not solely determined by the internal structure of the nu- 
cleon. Then, they also contain dynamical effects due to relativistic boosts and 
consequently the physical interpretation of the form factors becomes complicated. 
Recently, Kellv(|88|) has extracted spatial nucleon densities from the available form 
factor data. He selected a model for the Lorentz contraction of the Breit frame 
in which the asymptotic behavior of the form factors conformed to perturbative 
quantum chromo-dynamics (pQCD) scaling at large Q^-values and expanded the 
densities in a complete set of radial basis functions, with constraints at large 
radii. The neutron and proton magnetization densities are found to be quite 
similar, narrower than the proton charge density. He reports a neutron charge 
density with a positive core surrounded by a negative surface charge, peaking at 
just below 1 fm, which he attributes to a negative pion cloud. Alternatively, he 
extracts the radial distributions of the u and d quarks which both show a sec- 
ondary lobe which he interprets as an indication of an orbital angular momentum 
(0AM) component in the quark distributions. Friedrich & Walcher H89|) observe 
as a feature common to all EMFFs a bump/dip at Q ~ 0.5 GeV with a width of 
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~ 0.2 GeV. A fit to all four EMFFs was performed, assuming a dipole behaviour 
for the form factors of the constituent quarks and an Z = 1 harmonic oscillator 
behaviour for that of the pion cloud. They then transformed their results to co- 
ordinate space, neglecting the Lorentz boost, where they find that the pion cloud 
peaks at a radius of ~ 1.3 fm, slightly larger than Kelly did, close to the Compton 
wavelength of the pion. Hammer et al.(|9n)l argue from general principles that the 
pion cloud should peak much more inside the nucleon, at 0.3 fm. However, 
they assign the full NN2tt continuum to the pion cloud which includes different 
contributions than just the one-pion loop that Kelly (and Friedrich & Walcher) 
assign to the pion cloud. The structure at ~ 0.5 GeV, common to all EMFFs, is 
at such a small Q^-value that its transformation to coordinate space should be 
straightforward. 

2.4 Model Calculations 

The recent production of very accurate EMFF data, especially the surprising 
data from polarization transfer, has prompted the theoretical community 
to intensify their investigation of nucleon structure. Space limitations compel 
us to focus on only a few highlights. The interested reader is encouraged to 
read the original publications; the review by Thomas & Weise (|91() is an excellent 
introduction. 

The U-, d- and s-quarks are the main building blocks of the nucleon in the kine- 
matic domain relevant to this review. Its basic structure involves the three lightest 
vector mesons (p, lo and </>) which have the same quantum numbers as the pho- 
ton. Consequently, one should expect these vector mesons to play an important 
role in the interaction of the photon with a nucleon. The first EMFF models 
were based on this principle, called vector meson dominance (VMD), in which 
one assumes that the virtual photon - after becoming a quark-antiquark pair - 
couples to the nucleon as a vector meson. The EMFFs can then be expressed 
in terms of coupling strengths between the virtual photon and the vector meson 
and between the vector meson and the nucleon, summing over all possible vector 
mesons. In the scattering amplitude a bare-nucleon form factor is multiplied by 
the amplitude of the photon interaction with the vector meson. With this model 
lachello et al. (|92j) predicted a linear drop of the proton form factor ratio, similar 
to that measured by polarization transfer, more than 20 years before the data be- 
came available. Gari & Kriimpelmann i93| extended the VMD model to conform 
with pQCD scaling at large Q^-values. The VMD picture is not complete, as 
becomes obvious from the fact that the Pauli isovector form factor F2 is much 
larger than the isoscalar one F2 ■ An improved description requires the inclu- 
sion of the isovector vrvr channel through dispersion relations(|2^^2l))- By adding 
more parameters, such as the width of the p- meson and the masses of heavier 
vector mesons (|94() . the VMD models succeeded in describing new EMFF data as 
they became available, but with little predictive power. Figure |H1 confirms that 
Lomon's calculations provide an excellent description of all EMFF data. Bijker 
& Iachello ()97() have extended the original calculations by also including a meson- 
cloud contribution in F2, but still taking only two isoscalar and one isovector 
poles into account. The intrinsic structure of the nucleon is estimated to have 
an rms radius of ~ 0.34 fm. These new calculations are in good agreement with 
all EMFF data, except for at low Q^-values. The most recent dispersion- 
theoretical analysis (jTSJ, using four isoscalar and three isovector mesons, results 
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in an excellent description of and G^j, but only reasonably describes and 
Gg. Subsequent studies (j98jl have further developed this combined approach to 
include chiral perturbation theory. However, such models can only be used at 
smah Q^.values, < 0.4 GeV^. 




02460246 
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Figure 9: Comparison of various calculations with available EMFF data, indi- 
cated by the same symbols as in Figures |31 El El a-nd [Tj For only polarization- 
transfer data are shown. Not shown are the data for G'^ of References (|721I66|) 
and the data for Gl^ of References (•p.ip,'Sl,59,,56;). For G% the results of Schi- 
avilla &: Sick(j75|) have been added. The calculations shown are from References 
H97II94II78II11UII111II121|) . Where applicable, the calculations have been normalized 
to the calculated values of ^p,n- 

Many recent theoretical studies of the EMFFs have applied various forms of a 
relativistic constituent quark model (RCQM). Nucleons are assumed to be com- 
posed of three constituent quarks, which are quasi-particles where all degrees of 
freedom associated with the gluons and qq pairs are parametrized by an effective 
mass. Because the momentum transfer can be several times the nucleon mass, the 
constituent quarks require a relativistic quantum-mechanical treatment. Three 
possibilities exist for such a treatment: the instant form, where the interaction is 
present in the time component of the four-momentum and in the Lorentz boost; 
the point form, where all components of the four-momentum operator depend on 
the interaction; and the light-front form, where the interaction appears in one 
component of the four- momentum and in the transverse rotations. In each of 
these forms the Poincare invariance can be broken in the number of constituents 
(by the creation of qq pairs) or by the use of approximate current operators. 
Although most of these calculations correctly describe the EMFF behaviour at 
large Q^-values, effective degrees of freedom, such as a pion cloud and/or a finite 
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size of the constituent quarks, are introduced to correctly describe the behaviour 
at lower Q^-values. 

Miller (|99|1 uses an extension of the cloudy bag model HlUU|) . three relativisti- 
cally moving (in light-front kinematics) constituent quarks, surrounded by a 
pion cloud. He chose a spatial wave function, as derived by Schlumpf , 
whose parameters (and those of the pion cloud) are chosen to describe the mag- 
netic moments, the neutron charge radius, and the EMFF behavior at large 
Q^-values. Cardarelli & Simula(l in2|l also use light-front kinematics, but they 
calculate the nucleon wave function by solving the three-quark Hamiltonian in 
the Isgur-Capstick one-gluon-exchange potential. In order to get good agreement 
with the EMFF data they introduce a finite size of the constituent quarks in 
agreement HlU3|l with recent DIS data. The results of Wagenbrunn et al. (|lU4() 
are calculated in a covariant manner in the point-form spectator approximation 
(PFSA). In addition to a linear confinement, the quark-quark interaction is based 
on Goldstone-boson exchange dynamics. The PFSA current is effectively a three- 
body operator (in the case of the nucleon as a three-quark system) because of its 
relativistic nature. It is still incomplete but it leads to surprisingly good results 
for the electric radii and magnetic moments of the other light and strange baryon 
ground states beyond the nucleon. Although Desplanques and Theussl (;105)) have 
criticized the use of the point form in its introduction of two-body currents in the 
form of a neutral boson exchange, Coester and Riska ()lU6|) obtain a reasonable 
representation of empirical form factors in this frame. Giannini et al.fTWj) have 
explicitly introduced a three-quark interaction in the form of a gluon-gluon inter- 
action in a hypercentral model, which successfully describes various static baryon 
properties. Relativistic effects are included by boosting the three quark states to 
the Breit frame and by introducing a relativistic quark current. All previously 
described RCQM calculations used a non-relativistic treatment of the quark dy- 
namics, supplemented by a relativistic calculation of the electromagnetic current 
matrix elements. Merten et al. (|108p have solved the Bethe-Salpeter equation with 
instantaneous forces, inherently respecting relativistic covariance. In addition to 
a linear confinement potential, they used an effective flavor-dependent two-body 
interaction. For static properties this approach yields results (|in9)l similar to those 
obtained by Wagenbrunn et al. (|in4|) . The results of these five calculations are 
compared to the EMFF data in Figure 1101 The calculations of Miller do well 
for all EMFFs, except for at low Q^-values. Those of Cardarelli & Simula, 
Giannini et al. and Wagenbrunn et al. are in reasonable agreement with the data, 
except for that of Wagenbrunn et al. for G^j, while the results of Merten et al. 
provide the poorest description of the data. 

Before the Jefferson Lab polarization transfer data on G^/G^^ became available 
Holzwarth ( 110) predicted a linear drop in a chiral soliton model. In such a model 
the quarks are bound in a nucleon by their interaction with chiral fields. In the 
bare version quarks are eliminated and the nucleon becomes a skyrmion with a 
spatial extension, but the Skyrme model provided an inadequate description of 
the EMFF data. Holzwarth's extension introduced one vector-meson propagator 
for both isospin channnels in the Lagrangian and a relativistic boost to the Breit 
frame. His later calculations used separate isovector and isoscalar vector-meson 
form factors. He obtained excellent agreement for the proton data, but only 
a reasonable description of the neutron data. Kim et al. lTIII used an SU(3) 
Nambu-Jona-Lasinio Lagrangian, an effective theory that incorporates sponta- 
neous chiral symmetry breaking. This procedure is comparable to the inclusion 
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Figure 10: Comparison of various RCQM calculations with available EMFF data, 
similarl to the comparison in Figure IHl The calculations shown are from Refer- 
ences (piirn^lWllTnillTnHI) . Mlller (q-only) denotes a calculation by Miller (PI) 
in which the pion cloud has been suppressed. 

of vector mesons into the Skyrme model, but it involves many fewer free param- 
eters (which are fitted to the masses and decay constants of pions and kaons). 
The calculations are limited to Q^< 1 GeV'^ because the model is restricted to 
Goldstone bosons and because higher-order terms, such as recoil corrections, are 
neglected. A constituent quark mass of 420 MeV provided a reasonable descrip- 
tion of the EMFF data (Figure IHl). 

In the asymptotically free limit, QCD can be solved perturbatively, providing 
predictions for the EMFF behavior at large Q^-values. Brodsky & Farrar (|112j) 
derived a scaling law for the Pauli and Dirac form factors based on a dimen- 
sional analysis, that entailed counting propagators and the number of scattered 
constituents: 

Fioc(Q2)-2^ F2a(Q2)-3, F2/F10CQ-2 (16) 

Brodsy and LeDage H113|l later reached the same asymptotic behavior based on a 
more detailed theory that assumed factorization and hadron helicity conservation. 
The recent polarization transfer data clearly do not follow this pQCD prediction 
(which the Rosenbluth data unfortunately do). Miller pi4p was the first to ob- 
serve that imposing Poincare invariance removes the pQCD condition that the 
transverse momentum must be zero, and introduces a quark 0AM component in 
the wavefunction of the proton, thus violating hadron helicity conservation. His 
model predicts a 1/Q behaviour for the ratio of the Dirac and Pauli form factors 
at intermediate Q^-values, in excellent agreement with the polarization transfer 
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data for > 3 GeV^. Iachello (|77() and others have pointed out that this 1/Q 
behaviour is accidental and only valid in an intermediate Q^-region. Ralston ()115() 
has generalized this issue to conclude that the Q^-behavior of the Jefferson Lab 
data signals substantial quark 0AM in the proton. Recently, Brodsky et al. (|116j) 
and Belitsky et al.(117) have independently revisited the pQCD domain. Belitsky 
et al. derive the following large Q^-behavior: 

where A is a soft scale related to the size of the nucleon. Even though the Jefferson 
Lab data follow this behavior fFigure[TT|). Belitsky et al. warn that this could very 
well be precocious, since pQCD is not expected to be valid at such low Q^-values. 
Brodsky et al.(116) argue that a nonzero 0AM wave function should contribute 
to both Fi and F2 and that thus Q'^F2/Fi should still be asymptotically constant. 




Figure 11: The ratio {Q"^ F2 / Fi) / {Q"^ / A"^) as a function of foj. ^^e 
polarization-transfer data and the calculations of References (|97[l^[TTUlllU7j) . 
The same ratio, scaled by a factor -1/15, is shown for the neutron with open 
symbols. For A a value of 300 MeV has been used. 

Once enough data have been collected on generalized parton distributions, it 
will become possible to construct a three-dimensional picture of the nucleons, 
with the three dimensions being the two transverse spatial coordinates and the 
longitudinal momentum. Miller(j[llB| has further investigated the information 
that can be extracted from form-factor data by themselves. His colorful images 
of the proton should be interpreted as three-dimensional pictures of the proton 
as a function of the momentum of the quark, probed by the virtual photon, and 
for different orientations of the spin of that quark relative to that of the proton. 
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Ji H119|l has derived similar images from generalized parton distributions using 
Wigner correlation functions for the quark and gluon distributions. 
However, all theories described until now are at least to some extent effective 
(or parametrizations) . They use models constructed to focus on certain selected 
aspects of QCD. Only lattice gauge theory can provide a truly ab initio calcu- 
lation, but accurate lattice QCD results for the EMFFs are still several years 
away. One of the most advanced lattice calculations of EMFFs has been per- 
formed by the QCDSF collaboration!^^. The technical state of the art limits 
these calculations to the quenched approximation (in which sea-quark contribu- 
tions are neglected), to a box size of 1.6 fm and to a pion mass of 650 MeV. 
Ashley et al.fl21) have extrapolated the results of these calculations to the chiral 
limit, using chiral coefficients appropriate to full QCD. The agreement with the 
data (Figure ^ is poorer than that of any of the other calculations, a clear indi- 
cation of the technology developments required before lattice QCD calculations 
can provide a stringent test of experimental EMFF data. 



3 Generalized Polarizabilities of the Nucleon 



3.1 Introduction 

The electric and magnetic polarizabilities of an object describe how its internal 
structure responds to external electric and magnetic fields. In the weak field 
limit, an external electric field E induces in a finite system an electric dipole 
moment p proportional to the applied field: 

p = ueB (18) 

(Heaviside-Lorentz units). This proportionality defines the electric polarizability 
oe- The induced dipole moment p is measurable by the effects of the long range 
dipole electric field it produces. Similarly, a weak external magnetic field H 
induces a magnetic dipole moment 

A/2 = /?mH, (19) 

which defines the magnetic polarizability Pm- The induced moment A/2 is a 
change in the static moment /io of a non spin zero system. 

It is instructive to contrast a few basic examples of the polarizabilities. For the 
hydrogen atom, the electric polarizability is of the same order of magnitude as 
the hydrogen atom volume, whereas the proton's polarizability is roughly 10"'^ 
of its volume. Thus the proton is much stiffer than the hydrogen atom. In the 
Schrodinger equation for the hydrogen atom, the electric polarizability is ()122|) : 



atom _ 27 

^ Svr 



4 3 
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where ao = l/{meaQED) is the Bohr radius. For a particle of mass m and 
charge e bound by a simple harmonic oscillator of length constant b, the electric 
polarizability is 
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where Ac = l/m is the Compton wavelength of the particle. The Hydrogen 
atom result is compatible with a non-relativistic harmonic oscillator model with 
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6 ~ ao and therefore b/Xc ^ 1- On the other hand, a harmonic oscillator 
model of the small value of the proton polarizability requires b/Xc ~ 1 for the 
quarks or pions of the proton substructure. Thus the order of magnitude of 
the nucleon polarizabilities is evidence for the explicitly relativistic structure of 
the nucleon. This illustrates how the polarizabilities and their dependent 
generalizations reveal details of the nucleon dynamics that go beyond the charge 
and magnetization distributions of the form factors. 



3.2 Proton Polarizabilities 



The low-energy limit of the Compton amplitude is determined by the nucleon's 
charge and magnetic moment. As the energy of the probe increases, the effect 
of the excitation spectrum of the nucleon can be summarized by a set of electro- 
magnetic polarizabilities. Today there is an extensive experimental program to 
measure the proton polarizabilities in real and virtual Compton scattering and 
the neutron polarizabilities in elastic and quasi-elastic Compton scattering on the 
deuteron. 

The scattering amplitude for elastic real photon scattering from a nucleon can 
be described by six complex amplitudes, each of which multiplies a linearly in- 
dependent algebraic structure of the polarization and kinematic variables (|123)1 . 
Below pion production threshold these amplitudes are real, to lowest order in 
ctQED = e^/(47r). 

To illustrate the low energy structure of Compton scattering, consider first the 
forward Compton amplitude on a nucleon of charge Ae and anomalous magnetic 
moment k. The scattering amplitude has the following form ()1241ll25|l : 
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Here u is the laboratory photon energy, e* and e ' the initial and final photon polar- 
ization vectors, and a the nucleon spin operator. The forward spin polarizability 
7o will be discussed later. 

To order a;^, the spin-averaged Compton cross section is given by (|124|ll25l> : 
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da^ /da^y is the exact (Born) cross section for a nucleon, given by Powell H126|) . 
It differs from the Klein-Nishina formula fiZT)) in the inclusion of the anomalous 
magnetic moment, k. The linear dependence of the cross section on the polariz- 
abilities oe and (3m at order u)"^ results from the interference of the ©(u;^) ampli- 
tude with the order unity Thomson amplitude. For Compton scattering from a 
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free neutron, since A = 0, this interference is absent and the polarizabihties enter 
the cross section at order u;^. 

Each of the six independent terms in the Compton amphtude can be constrained 
by dispersion relations (DR). For the forward Compton amphtude, these relations 
take a particularly simple form, since the optical theorem connects the imaginary 
part of the forward scattering amplitude to the photo-absorption cross section. 
For /2, the dispersion relation yields the Gerasimov-Drell-Hearn (GDH) sum rule 



aQED 2 ^ 



Th 



I du! 

0"3/2(^*^) - (^l/2{^)\ — (25) 



For /i, the subtracted dispersion relation yields the Baldin sum rule ()13Uj) 

1 POO ~ 

aE + pj^ = — / -^^dw = [14.20 ± 0.50] • 10"^ fm^ (1970), (26) 

27^ JTh W 

= [13.69 ±0.14] • lO-'^fm^ (1998). (27) 

The l/w^ behavior of the integrand of Equation 1261 gives a rapid convergence of 
the integral. Already in 1970, Damashek & Oilman ()131|1 were able to obtain 
the value quoted in Equation 1261 The intervening 30 years of experiments have 
resulted in the improved analysis of Babusci et al. (132) in Equation 1271 This 
convergence is to be contrasted with the much more difficult problem of the ex- 
perimental evaluation of the GDH sum rule (Equationl25|). since the experimental 
integration saturates only if a model is used for the (unmeasured) high energy 
behavior of the integrand. In contrast with Equation 1261 the dispersion relation 
for aE — Pm is not convergent. Thus, this combination of polarizabihties must be 
determined directly from low energy Compton scattering. In an effective theory, 
the lowest order interaction of the external electromagnetic field with the inter- 
nal structure of the proton (beyond the charge and magnetic moment) is entirely 
described by the static electric and magnetic polarizabihties: 



n 



(2) 
Int 



'i'KUEF? + A-k(3m^^ ■ (28) 



L'vov showed that the 0{uj'^) terms in the scattering amplitude are the same 
polarizabihties (|133|) . 

Starting with Gol'danski (1134() . several generations of tagged and untagged bremsstrahlung 
experiments have tackled the measurement of low energy Compton scattering 
from the proton (I1351I1361[T?71[T!M[T!P) : see (|Tin|l for a historical review. The 
bremsstrahlung photon beam is obtained by passing the primary electron beam 
through a high Z radiator (such as Cu or W), a few percent of a radiation 
length in thickness. The electron beam is deflected by a magnet so that only 
the photons reach the hydrogen production target. In the tagged experiments, 
the bremsstrahlung electrons are detected in coincidence with the scattered pho- 
ton, thereby "tagging" the incident photon energy of the Compton event. In the 
untagged experiments, the incident photon energy is determined event-by-event 
from the constrained kinematics of the elastic scattering process. 
The most recent and precise Compton results are from the MAMI accelerator 
H141|) . In order to extract the polarizabihties from Compton scattering data at 
and above pion threshold, the higher order terms in the full scattering ampli- 
tude are constrained by dispersion relations, derived in particular by L'vov, and 
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collaborators ()142lll4t^ . The dispersion relations are evaluated by pion photo- 
production multipole analysis along with a Regge theory extrapolation for the 
asymptotic part. Drechsel et al. 1)14411145'}) also discuss subtracted dispersion 
relations for two of the six amplitudes. 

A global analysis of the Compton data yielded the following values for the polar- 
izabilities (|141j) : 

aE = [12.1 ± 0.3{stat.) =F OA{syst.) ± 0.3{model)] ■ W'^ fm^ 

Pm = [1-6 ± OAistat.) ± OAisyst.) ± OA{model)] ■ 10"^ fm^. (29) 

The systematic errors are anti-correlated from the Baldin sum rule constraint 
(Equation |27| . The model dependent errors are estimates of the uncertainties 
from the dispersion relations. 

The 0{uj'^) spin-dependent and 0{lo^) terms in the Compton scattering amplitude 
can also be connected to polarizabilities, as discussed in detail by Babusci et al. 
p25jl . In the presence of time- or spacedependent electromagnetic fields, the 
internal structure of the nucleon will polarize in response to the form of the 
external field. The most general effective interaction (to order u!^) between the 
external field and the internal nucleon structure is given by (|125)1 : 
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(31) 

where Eij = i^iEj + V-,£'j)/2 and similarly for Hij. Because of the extra space 

and time derivatives in Equation |^ the polarizabilities defined by and 
enter the spin-dependent and spin-independent Compton scattering amplitudes 
to order and w^, respectively. The polarizabilities -)e2 and 7m2 measure the 
spin dependent quadrupole strength in the nucleon spectrum. DR estimates 
and theoretical calculations of these higher order polarizabilities are presented in 
References ((T^ and (fT2'2| . 

Of particular experimental interest are the forward 70 and backward 7,^ spin 
polarizabilities: 

70 = -lEl - IMI - 1E2 - 1M2 (32) 

= -^y ^ [^^3/2('^)-^^i/2('^)] dnniiiizi) (33) 

l-n = -JEI + IMI + 1E2 - 1M2 (34) 

The forward spin polar izability 70 is the 0(u;^) term in the forward scattering 
amplitude of Equation [231 Similarly, '^^^ is the 0(a;^) spin-dependent term in the 
amplitude for Compton scattering in the backward direction (Q^^ = vr). 
Sandorfi et al. extracted the forward spin polar izability from a dispersion analysis 
of photoproduction multipoles (PH|l . The MAMI (|T¥?|l and ELSA ((T^ GDH 
experiments, measured the following contributions (units 10~^ fm^), respectively, 
70 = -1.87 ± 0.08 ± 0.10 for 0.2 GeV < u; < 0.8 GeV and 70 = -0.027 ± 0.002 ± 
0.001 for 0.8 GeV < a; < 1.82 GeV. From this, along with the multipole estimate 
70 = +0.90 for u < 0.2 GeV (I14'^1I151)) one obtains: 

70 = [-1.02 ± 0.08(stat) ± OAO{syst)] ■ lO'^fm^ (35) 
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There are several measurements of the backward spin polarizabihty 77r. The over- 
all error envelopes in the analysis of the low energy Compton data (Equation 129(1 
increase slightly if is included as a fitting parameter, with the result ()141|) : 

7^ = [-36.1 ± 2.1{stat.) =F 0.4(syst.) ± 0.8{model)] ■ 10"^ fm^ (36) 

The LEGS group obtained the result (|1 52111 r>:^ 

= [-27.1 ± 2.3{stat. + syst) ± 2.2{model)] ■ fm^ (37) 

from back-scattering jp 'jp and jp — > jp data in the A-resonance region. 
However, new Compton cross section measurements in the A-resonance confirm 
the larger value (,154tll55)) 

7^ = [-37.1 ± 0.6{stat. + syst) ± 3.0{model)] ■ 10"^ fm^. (38) 

This result was again confirmed with a second apparatus at MAMI ()156|) . When 
one compares the results of Equations OH and |3H1 the high energy data yield an 
improved statistical precision at the expense of a greater model uncertainty. The 
model uncertainties include distinct analysis with the MAID p57j) or SAID p59j) 
multipoles, and variations within each multipole parameterization. 
The 7771^^ triangle anomaly dominates the backward spin polarizabihty via the 
vr'^ t-channel exchange(16lT): 

7^^") = -45.0 ± 1.6 • 10-^ fm^. (39) 

This contribution is absent from 79. The larger experimental value of |7^| is con- 
sistent with DR calculations (jl6npi61|) . In addition to the vr'^ anomaly, L'vov and 
Nathan { 161^) compute contributions of -1-7.3, —0.3, and —1.6 (units of 10~^ fm'') 
from vrA^ s-channel, vrvrA^ s-channel, and rj, ry' t-channel contributions, respec- 
tively. Their total DR prediction is 7,^ = (-39.5 ±2.4) • 10^ fm^. This calculation 
also illustrates that in contrast with aE — Pm, for which Regge phenomenology 
suggests the dispersion relations do not converge, or oe + Pm, for which the vrvrA^ 
channel contributes 15% to the Baldin sum rule, the dispersion relations for the 
higher order polarizabilities converge rapidly. 



3.3 Theoretical Perspective 

In Heavy Baryon Chiral Perturbation Theory (HBChPT) to 0{p^), the polariz- 
abilities have a very simple form (|122lll6li)) : 
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where gA_ ~ 1.266 is the nucleon axial coupling constant and F^^ ~ 92.4 MeV is the 
pion decay constant. The HBChPT results are in remarkable agreement with the 
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proton data. Even for 7q, if the tt" anomaly term 12^o/gA in 7,r is considered to 
set the scale for each of the four terms in Equation 1331 the disagreement between 
theory and experiment is small. 



Bernard et al. (^HIE^ obtained the following values for oe and /3m of the 
nucleon to 0{p'^) in HBChPT: 



The theoretical error bars come from the phenomenological constants in the ef- 
fective Lagrangian, especially the terms governing the A-Resonance contribution. 
Hemmert et al. introduced the A as an explicit degree of freedom in the 

the small-scale expansion (SSE) in which the mass difference M^ — Mjy is treated 
as an expansion parameter to 0{p^) along with ttItt and the momentum p. They 
obtain a large (4 — 7 • 10~^ fm^) additional positive contribution to a^; from Avr 
loops and to Pm from the A-pole (see below, Equation EHl) • These effects are ex- 
pected to be canceled in higher orders, as suggested in Reference (,163)1 . For the 
spin polarizabilities 7_bi,mi,E2,M25 the 0{p'^) ChPT calculations do not introduce 
new free parameters (|166|ll67lll68r i69, 170^. However, the values for 70 and 7,r 
change by approximately the magnitude of 70, relative to the 0{p^) results. 
These theoretical results leave somewhat in question the convergence of the chiral 
expansions for the polarizabilities. However, combined with the data, they con- 
firm basic expectations. The magnetic polarizability has a strong cancellation 
between paramagnetism (from the constituent quark spin flip, or equivalently 
— > N* terms) and diamagnetism driven by the pion cloud. For the electric 
polarizability, both degrees of freedom contribute with the same sign. It remains 
for the study of generalized polarizabilities to see if these contributions have a 
different spatial structure. 

Instead of attempting to predict the polarizabilities, Beane et al. H171|l used the 
electric and magnetic polarizabilities as the only two free parameters in a 0{p^) 
chiral perturbation expansion of the Compton amplitude H172|) . After refitting 
the data below 200 MeV, they obtain 



Magnetic polarizabilities for hadrons were calculated in Lattice QCD by Zhou 
et al. (jTTB). for values of 0.3 < < 1 GeV^. They added a static magnetic 
field to the lattice, and extracted the polarizability from the quadratic depen- 
dence of the ground state mass on the external field (Eauation l28|) . Substantially 
greater computational resources are needed to obtain results closer to the physical 
pion mass, but at the scale = 0.3 the results are = ± 1 and (3^^^ = 11 it 1. 
Christensen et al. (| 174(1 used the same technique to obtain electric polarizabil- 
ities of neutral hadrons. The result for a"^ shows considerably larger numerical 
uncertainty than the result. 

3.4 Neutron Polarizabilities 



= (10.5 ±2.0) •lO^'^fm^ 
Oe = (13.4 ± 1.5) • lO-'^fm^ 



Plj = (3.5 ± 3.6) • 10^^ fm^ 

= (7.8 ± 3.6) • 10"^ fm^ (44) 




(Chiral Fit) 



(45) 



consistent with the DR analysis of the same data, but with slightly larger uncer- 
tainties. 



Wissmann et al. reviewed the experimental techniques for measuring the neutron 
polarizabilities (|175|1 . Two recent analysis of the total photo-absorption data on 
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the deuteron and (7, vr) threshold multipoles give the following values for the 
Baldin sum rule of the neutron polarizabilities: 

al + Pli = 14.40 ± 0.66 • 10"^ fm^ (46) 
al + p^i^ = 15.2 ±0.5- 10"^ fm^ ^EB- (47) 

Schmiedmayer et al. (|176() and Koester et al. p77|) extracted the neutron electric 
polarizability from the energy dependence of low energy neutron-nucleus scatter- 
ing. Their results are: 

an = [12.6 ± 1.5 ± 2.0] • 10"^ fm^ (fTTHj) (48) 
an = [0.6 ± 5] • 10"^ fm^ ((nZ|). (49) 

Koester et al. (IT77') and Enik et al. (|178)1 suggest the uncertainty in both 
measurements should be ±5 • 10~^ fm^. The neutron polarizabilities were also 
extracted from D^j, and D(7, ^n)p measurements at MAMI-A (|179pi80|) the 
Saskatoon Accelerator Laboratory (SAL) H1811I182|) . Lund MAX-Lab (fTHH|l and 
MAMI (ITHUl . 

The most precise quasi-elastic D('y,n)p data were obtained at MAMI, with the 
result (UHIIHSI) 

a%-(3l[ = 9.8 ± 3.6(stat) "^J'l (syst) ±2.2(mode0 • 10"^ fm^. (50) 

The Mainz experiment also extracted the polarizabilities of the bound proton via 
the D{'y,^p)n reaction, with the result (|186j) : 

[a^E - ^m]^"''""^ = 9.l±1.7{stat + sys)±1.2{model)-W-^im^. (51) 

Although they did not extract a value for 7,^ for a bound proton, the authors note 
that their analysis is consistent with the free value of = —37.6 and inconsistent 
with = —27.1 (units of 10~^ fm^). The agreement with the free proton values 
for aE — Pm and 'Jt^ confirms the basic validity of the theoretical framework for 
the extraction of nucleon polarizabilities from quasi-free Compton scattering on 
the deuteron (|187j) . 

Hornidge et al. (|181|) and Lundin et al. fl83|) measured the coherent Compton 
scattering on the deuteron: 'jD — > 'jD. The data from SAL are at the highest 
energy (94 MeV) and therefore have the greatest sensitivity to the polarizabilities, 
but they are integrated over a 20 MeV energy bin, compared to the 10 MeV bins of 
the Lund data. From a global analysis using the A^iV-potential model formalism 
of Levchuk & L'vov (188), Lundin et al. extract the isoscalar polarizabilities 

a^ + PZ = 16.7 ± 1.6- 10"^ fm^ 

a^-P^ = 4.8 ± 2.0 • 10"^ fm^ (52) 

The theoretical formalism for deuteron Compton scattering is also discussed in 
Reference (|189j) . As already described for the proton Compton case, Beane et al. 
reanalyzed the data in a chiral expansion of the scattering amplitude, and ob- 
tained (|171j) 



10-^ fm^ 
lO'^fml 



(53) 
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The error bars result from both model uncertainties and inconsistencies between 
the data sets, with the chiral analysis more consistent with the Lund data alone. 
Measuring neutron polarizabilities to the same precision as the proton remains a 
formidable challenge. 

3.5 Virtual Compton Scattering 




VCS Born Non-Born 



Figure 12: Kinematics and scattering amplitude for the ep — > ep'j reaction. The 
incident and scattered electron 4-momentum vectors are and k^' , respectively. 
The initial and final proton 4-momentum vectors are p and p' , respectively. The 
final photon and the VCS virtual photon 4-momentum vectors are q' and = 
{k^ — k^'), respectively. The Mandelstam invariants are s = {q+p)"^, t = {q — q')'^, 
and u = {p — q')"^. 

Virtual Compton Scattering (VCS) can be measured in the ep — > ep^ reaction. 
In this case, the Compton amplitude interferes with the Bethe-Heitler (electron 
radiation) amplitude, as illustrated in Figure 1121 which also defines the kine- 
matic variables. Guichon et al. constructed a gauge invariant separation of the 
ep — > ep7 amplitude into the Bethe-Heitler (BH), Born and Non-Born (NB) 
terms flQOl) . The BH and Born terms are the amplitudes for electron and proton 
bremsstrahlung, including only the on-shell proton form factors -Fi,2- They then 
expanded the amplitude in powers of g', the final photon energy in the proton- 
photon center-of-mass frame. The leading term in the expansion of the BH and 
Born amplitudes is 0{l/q'), arising from the electron and proton propagators, 
respectively. The leading term in the Non-Born amplitude is 0{q'). This term is 
defined by 6 generalized polarizabilities, representing the independent multipoles 
coupling the initial virtual photon with a final El or Ml photon H19U|) . The 
generalized polarizabilities are functions of Q^, the invariant momentum transfer 
squared in the ^ limit. The generalized polarizabilities describe the spatial 
variation of the polarization response of the proton, as described explicitly in 
Reference (fTM|l . 

Metz & Drechsel applied crossing and charge-conjugation symmetry in the linear- 
cr model, and obtained four constraints among the ten low energy VCS multipoles 
(|192U19^ . Later, Drechsel et al. showed that these constraints are general. 
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establishing that there are six independent generahzed polarizabihties ()194II195|) . 
This is a nice apphcation of model building: an approximate model incorporating 
chiral symmetry and exact relativistic dynamics led to a deeper understanding 
of the fundamental dynamics. 

The generalized polarizabihties are labeled by P^p ^ '^^^^{Q"^), where L' and L 
denote final and initial multipolar ity, respectively; p' and p indicate the polar- 
ization of the final and initial photon, respectively, which may be Coulomb (C), 
magnetic (M), or electric (E); and 5 = 0, 1 for a scalar or vector operator in nu- 
cleon spin space. Siegert relations connect the electric and Coulomb multipoles 
096 ). The independent set of generalized polarizabihties that enter the cross 
section to lowest order (beyond the BH-|-Born terms), and their ^ limits, 
are Q96j): 



p(Cl,Cl)l(g2) ^ p(Ml,A^l)l(g2) ^ q 



aQEDP'^''''''''>\Q') - -(2/3)3/2 



7-B2 



aQEDP^'''''''^'m - -{V2/3)^M2 (54) 
In an unpolarized VCS experiment the cross section (to order q'^) has the form: 

da = da'''^+''°''' + VLL[PLL-PTT/e]+VLTPLT, (55) 

PlL = -V24MG^(g2)p(Cl.^l)0(Q2)^ J^GP(g2)^^(g2) (5g) 

OIQED 

Ptt = 6M(l + f)G^(Q2)[p(Mi.Mi)i(Q2)^^~p(ci,Af2)i(Q2)] (57) 

= yiMTTTf [G^(Q2)p(*^^'^'^i)°(g2)-V6G^,(Q2)p(^^'^i)^(g2)" 

2^ . VTTf G^(g2)/3^(Q2) _ gpi^ (58) 



OiQED 



where vll and vtt are kinematic factors defined in Equations 97-100 of Refer- 
ence 1^^, e = 1/[1 -|- 2(q^/(5^) tan 0/2]iab is the virtual photon polarization (in 
the VCS amplitude), and f = Q"^ /{4M^). The generalized polarizabihties enter 
linearly in the cross section to this order, owing to the interference between the 
NB and BH-hBorn terms. In the ^ (RCS) limit, Pll ^ AMaE/aQED, 
Ptt 0, and Plt -'iM(5M/oiQED- 

Equation 1551 illustrates that an unpolarized VCS experiment (including no vari- 
ation in photon polarization) can measure only two linear combinations of the 
polarizabihties. The spin polarizabihties ^e2 and ^m2 enter the unpolarized (NB) 
VCS amplitude to lowest order even though they only enter the spin dependent 
RCS amplitude to order (Equation [231) • First, ^e2 and 7m2 enter the un- 
polarized VCS cross section because the virtual photons are linearly polarized, 
and this induces a polarization of the proton in the interfering BH-|-Born term. 
Second, in the TL^ terms of Equation |^ if the space derivatives are assigned to 
the initial virtual photon, these terms enter the VCS amplitude to the same order 
in q' as do the ordinary polarizabihties. 

Roche et al. 097| measured the VCS cross section on the proton below threshold 
at MAMI at fixed qcA/ = 600 MeV/c (Q^ = o.33 GeV^) and e = 0.62, where 
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Figure 13: MAMI ep — > ep^y data at fixed |q| = 0.6 GeV and five values of q' (botli 
variables in the photon-proton center-of-mass frame) H197|) . The differential cross 
sections are plotted as a function of the polar center-of-mass angle 9 between 
the virtual photon direction q = k*^ — k*^' and the outgoing photon q'. Values of 
6 < refer to kinematics with the azimuth of q' around q equal to 180°. The 
solid curves are the BH+Born calculations. The dashed curves include the low 
energy expansion of Equation 1551 with the two structure functions Pll — Ptt/^ 
and Plt fitted to the data. 



qcM is the VCS virtual photon 3-momentum in the proton-photon center-of- 
mass frame. The cross sections were extracted including radiative corrections 
calculated specifically for the full VCS process (fmsfl . Figure CHI displays the 
angular distributions of the ep — > ep^y cross section for five values of q'. The 
rapid rise in the cross sections for 6 > is the tail of the first of the two BH 
peaks, when the radiated photon is parallel to the incident (first peak) or scat- 
tered electron direction. The broad rise in the cross section for 6 < —90° has 
a strong contribution from the Born term, which is approximately a boosted 
Larmor dipole-radiation pattern from the proton. The cross section at low q' is 
consistent with the BH+Born cross section within the 2.5 % experimental uncer- 
tainty, including the uncertainty in the elastic proton form factors. The deviation 
of the data from the BH+Born cross section grows linearly with g', as expected 
for the polarizabilities. 

The generalized polarizabilities have been calculated in the Constituent-Quark 
model (190, 199, 200)) . a tree- level Lagrangian model H201|l . the Skyrme model 
d^nSl), the Linear cr-model (|1921I193|1 and to 0{p^) Chiral Perturbation Theory 
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Figure 14: Generalized polarizability structure functions extracted from jp — >■ 7p 
((Bll, and ep ep-f data of MAMI (197) and Jefferson Lab (: 208 ) . The Jefferson 
Lab points are the low energy expansion and dispersion-relation (DR) set lb 
points at = 1 and the DR set II point at = 2. The plots are Pll — Ptt/^ 
(left) and Plt (right). The dashed, solid, and dotted curves represent calculations 
in 0{p^) Heavy Baryon Chiral Perturbation Theory (HBChPT) (206j). Dashed 
curves, left and right are the total HBChPT predictions, with e = 0.62 (MAMI 
value) at left. The solid curves at left and right are the HBChPT contributions 
of and /3m, respectively. The dotted curves are the HBChPT contributions 
of the spin polarizabilities Ptt/^ (left) and (right). The dot-dashed 

curves are the DR predictions for the same spin polarizabilities (209). All data 
and curves are divided by the Brash parameterization of G^{Q'^) ()51)1 . 

(ChPT) H204ll205ll206|) . The generalized spin polarizabilities have also been cal- 
culated in 0(/) HBChPT (tWIl . 

Figure El displays the polarizability structure functions Pll — Ptt/^ and Plt 
extracted from the MAMI VCS experiment, along with the RCS results ()141() 
and the Jefferson Lab VCS results (208). The figure also shows the 0{p^) (one 
loop) HBChPT results H206|l . and the DR predictions of the spin polarizabilities 

Hemmert et al. calculated the generalized polarizabilities to 0{p^) (one loop) in 
HBChPT and to ©(e^) in the small-scale -expansion (SSE) with Ma — M taken 
as a third expansion parameter (as well as 771,7^ and p) (|206j) . The HBChPT 
calculation yields analytic expressions for the generalized polarizabilities, which 
are plotted in Figure ITU The agreement between the calculations and the data 
at low is striking. 

The HBChPT calculation of Pm{Q'^) has the dramatic feature of rising at low q. 
As expected in a naive picture, this results, from a partial cancellation between 
the diamagnetism of the pion cloud and the paramagnetism of the core. Larger 
values of q probe shorter distance scales, and are therefore dominated by the 
paramagnetism. Eventually the finite size of the proton imposes the decrease of 
both the para- and diamagnetic contributions. 

The strong cancellation between para- and diamagnetism is also emphasized by 
both the SSE O(e^) calculation (|206j) and the tree-level effective Lagrangian 
model ()201|) . The SSE calculation agrees with the HBChPT calculations for 




Figure 15: Generalized polarizabilities oe and Pm of the proton. The (model- 
dependent) data points are extracted from Figure El by subtracting the DR 
predictions of the spin polarizability contributions from the data (Eqs. I57l58j) . 
The heavy solid curves ending at = 0.5 GeV^ are the HBChPT predictions 
for a E and /3m (same as Figure ITU) (f2fl6| . See text for discussion of the DR 
curves. 



the variation of the Generalized Polarizabilities, and for the magnitude of the 
spin-polarizabilities. However, as discussed above in the RCS section p.^f) . at 
the photon point: 

a[SSE 0{e^)] - a[HBChPT 0{p^)] = 4.2 • 10"^ 

P[SSE 0{e^)] - l3[HBChPT 0{p^)] = 7.2 • 10"^ (59) 

The large value of /3m in the SSE calculation comes from the paramagnetism of 
the ^ A transition. It is expected that this will be canceled by 0{e^) dia- 
magnetic terms. Similarly, in the tree-level effective Lagrangian model, there is 
a strong cancellation between the — > A paramagnetism and the diamagnetic 
contribution from higher resonances H201|l . 

Pasquini et al. developed a DR formalism for the VCS amplitude up to the 
Ntttt threshold(209). The imaginary part of the VCS amplitude is expressed 
explicitly by unitarity in terms of the 7*A^ — > vrA^ — > MAID multipoles 
(|1571I158|) . The real part of the amplitude is expressed as a dispersive integral 
(as a function oi v = {s — u)/(4M) at fixed t and Q^) over the imaginary part, 
by the Cauchy theorem. If the dispersive integrals do not saturate at finite 
energy, then an (^--independent) asymptotic piece is added to the amplitude. 
This represents, equivalently, either a semicircular contour in the complex v- 
plane or the contribution of channels beyond vrA^. Of the 12 VCS amplitudes, 
FiiQ'^i'^it), the dispersive integrals converge in principle for all but Fi and F^, 
based on Regge phenomenology. 

The asymptotic contribution to F^ is obtained from i-channel vr'^-exchange, in 
accord with the calculation of the backward spin polarizability 7,^ in RCS. The 
asymptotic contribution to -Fi is obtained from i-channel cr-exchange, with a 
phenomenological 7*70" vertex A/3((5^) that must be fitted to the VCS data. 
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In this approximation, the non-Born (NB) contribution to the amphtude Fi is 
expressed as: 

Ff^{Q^u,t) = F[''{Q\u,t)+F r{Q',0 ,t) 

= Fi [Q ,u,t) + d , (60) 



where E = \/ + is the initial proton energy in the photon-proton center-of- 
mass frame in the g' — > hmit. The model dependence arises from the assumption 
that the asymptotic term F^^ is independent of v (at least below A^vrvr threshold). 
The magnetic polarizability is obtained from the (z^, f) = (0,— Q^) limit of the 
F^^ amplitude: 



MQ') = F^''{Q'MaQED^{E + M)/{2E). (61) 

The dispersive integral for F2 converges in principle, but in practice is poorly 
saturated by the MAID vriV multipoles (209). 

The asymptotic part of F2^ is the only contribution to [a + /3]((5^) that is not 
predicted by the DR calculations. In the absence of a multipole decomposition of 
the — > TTvrA^ amplitudes, the DR analysis approximates the low energy behav- 
ior of F2 in terms of the ttN multipoles plus a v- and f-independent asymptotic 
term: 

F^^{Q\v,t) = F2'^^(Q2,^,t) + F2^^^^0,0) 

Fr(Q^^,t) = F^(Q^^,t)-,P^ ' ,^[" + ^](^^).(62) 

^ ' ' ^ 2 w > , ; YF + M4M2(l + r) aQED 

The polarizability sum is f (|2U9() . Equation 29): 

[a£; + /?Ad(Q') = A[a + /3](Q2) 

-iM^aQED^I^^ [(1 + f)F2 + 2Fe + Fg - FuV"" {Q\ 0, -Q^) (63) 

In summary, the DR formalism of Reference H209() has a complete prediction of 
the VCS amplitude up to N-kit threshold, including all spin polarizabilities, in 
terms of just two unknown functions to be extracted from the data: A/3(Q2) and 
A[q-|-/?](Q2). At the real photon point, if one combines the DR calculations with 
the experimental results in Equation 1^ (|2U9() : 

f3l^{Q) = +9.1-10-^fm3 

A/3(0) = -7.5-10~^fm3 (64) 

[a + /?]^^(0) = 11.6- 10-^fm=^ 

A[a + /3](0) = +2.1-10-^fm3 (65) 

The vrA^ DR contribution is seen to be strongly paramagnetic, which arises natu- 
rally from the paramagnetic response of the constituent quarks which define the 
resonance spectrum. The asymptotic piece A/? is required to be strongly diamag- 
netic, which again is a natural result given the formal link with the pion cloud of 
the nucleon (via the t-channel a = [vrvrjo exchange). Finally, the asymptotic F2 
term A[a -|- /3](0) is only 15% of the total Baldin sum rule fEauationl27|). 
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The Jefferson Lab VCS collaboration analyzed ep — > ep7 data below pion thresh- 
old in terms of the low-energy expansion fEauation 155(1 and data up through the 
A-resonance in terms of the DR formalism. The results are shown in Figures ITU 
and 1151 Although the amplitudes Fi and F2 are the natural degrees of freedom 
in the DR formalism, the Jefferson Lab DR analysis follows (|2n9(l in estimating 
the asymptotic terms with two dipole form factors: 



AaE 



Aa(0) 



[1 + QV^l] 



2 



A,,, = (66) 

with the dipole parameters and A^ fitted to the data at each point. The 
dipole form is not essential to the analysis, since it is used only to describe the 
Q^-dependence within the acceptance of one spectrometer setting. 
In Figure El the DR predictions of the spin polarizabilities are used to extract 
o-EiQ"^) and Pm{Q'^) from the data. These data points are therefore subject to 
confirmation of the DR spin polarizabilities. In these two plots, we show the 
individual contributions to the total DR calculations. In the right-hand plot, the 
dot-dashed line is the contribution to of the vrA^ multipoles. The dashed line 
is the contribution of the phenomenological A/3 term, with A^ = 0.63 GeV. The 
thin solid line is the complete DR calculation for /3a/(Q^)- In the left-hand plot, 
the dot-dashed line is the contribution to aE of the nN multipoles (this has a 
small dia-electric contribution for > 0.1 GeV^). In accord with Equations 1611 
and 1631 we write the asymptotic contribution as AaE = A[a + (3] — A/3. The 
dashed curve is the contribution of the F^-^ term —A/3 to oe- Thus the pion 
cloud (i-channel cj-meson exchange) makes a diamagnetic contribution to /3m and 
a positive contribution to a^;. The dotted curve is the contribution of the F^-^ 
term, using a dipole form adjusted to the Jefferson Lab data: 

A[« + /3](g^) ^ + 



= 0.9 GeV (67) 

The dipole parameter A^ = 0.63 GeV is much smaller than the standard dipole 
fit to the nucleon form factors. This supports the interpretation that the pion 
cloud contribution has a much larger spatial size in the polarizabilities than in 
the form factors. The dipole form for A [a -|- /3] is unable to reproduce all of 
the data for ue- In particular, the DR curve falls well below the MAMI VCS 
point. The Sii(1535) and -Di3(1520) resonances have strong El couplings to 
the r]N and vrvrA^ channels, respectively, which are not included in the MAID 
analysis (|157|l209l) . Although the helicity amplitude for the -D13 resonance 
falls rapidly with Q^, the helicity amplitude rises by a factor of two from 
= to 1.0 GeV^ and the ^1/2 helicity amplitude of the 5ii resonance falls by 
less than 20% from Q2 = to 1.0 GeV^ (I15^[l21()[l^nil . Thus it is very plausible 
that the r]N and vrvriV channels contribute strongly to [a + /3] at short distance, 
as suggested by the data in Figure El 

The M.I.T.-Bates VCS (|212|) experiment used the Out-of-Plane-Spectrometer 
(OOPS) system to measure the full azimuthal distribution of recoil protons 
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Figure 16: Single and double spin asymmetries in the kinematics of the MAMI 
VCS experiment. Left: Beam helicity asymmetry in H(e, e'p)7. The four 
curves are the full DR predictions, with the asymptotic terms parameterized by 
(A«, A^) = (1.0,0.6) GeV (solid); (1.0,0.4) GeV (dashed); (1.0,0.7) GeV (dotted); 
and (1.4,0.4) GeV (dot-dashed). Right two panels: double-polarization observ- 
ables in H(e, e'p)'j, Z-axis parallel to q, X-axis in electron scattering plane. The 
dotted curves show the Bethe-Heitler-|-Born contribution; the solid curve shows 
the full DR calculation with (A^, A/3) = (1.0,0.6) GeV (same model as left panel); 
the dashed curves represent the HBChPT predictions H206|) . 

around the q direction at = 0.05 GeV^. This gives a strong sensitivity to 
the polarizabilities. This very low point will be valuable in testing chiral dy- 
namics and assessing the diverse length scales present in the proton polarization 
response. 

The HBChPT calculations in Figure suggest that the spin-polarizabilities make 
a large contribution to the unpolarized VCS observables, although the DR pre- 
dictions indicate otherwise. This emphasizes the need for double polarization 
measurements, that can directly measure the spin polarizabilities. These observ- 
ables were calculated by Vanderhaeghen (20^ , and are displayed in Figure 
The Mainz VCS collaboration has completed data taking on both single and dou- 
ble spin observables in H(e, e'p)7 and il{e,e'p)'y, in the kinematics of Figure IT^ 
The spin-polarizability term Ptt (Equation I55() can be isolated with a conven- 
tional Rosenbluth separation. This was discussed in the original Jefferson Lab 
VCS proposal ()213() . and will be feasible at MAMI with the energy upgrade. 
Laveissiere et al. show that higher measurements of the polarizabilities are 
feasible with Jefferson Lab at 6 GeV and above ff214,L 

3.6 Inclusive Electron Scattering and Forward Polarizabilities 

The inclusive electron scattering cross section on a nucleon has the form ()145() : 



da 



(68) 
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dT 



(TT + eai - hP^ J 2e{l - e) air - hP.Vl 



aQED 



rj^e/i lab 



K 1 



(Ttt (69) 
(70) 



In these expression, and k^' are the incident and scattered electron energies 
in the lab frame, and Pz and Px are the target nucleon polarizations parallel 
and perpendicular to q (Figure IT^ in the electron scattering plane. The virtual 
photon flux dT is evaluated in the Hand convention ()215|) : K = (s — M^)/(2M) = 
v — Q'^/{2M), with u = k'^ — k^' (in this sub-section) the electron energy loss. 
The partial cross sections are related to the usual DIS structure functions and 
the helicity cross sections (Jij2 and (T3/2 as follows. 



(JT 



0"L 



(Ttt 



<7l±l/2 



MK 
MK 
MK 

ATT^aQED 



MK ^ 

^TT^OQED ' 

MK 



l + l/Y) {M/u)F2{xBi,Q' 

gi{xBi,Q'^) - 7^5'2(a;Bj, Q^) 



-Fi{xBi,Q^] 
'71/2 - 0-3/2 



Fl{xBi,Q'^) T gi{xBi,Q'^) ±7^fi(2(xBj,Q^ 



(71) 
(72) 



with 7^ = Q'^/i^'^, and xbj = Q"^ /{2Mv) the Bjorken momentum-fraction variable 
(|216|) . cri±i/2 = 0-3/2 , o"i/2 are the virtual photo-absorption cross sections with 
total photon-|-proton helicity 3/2 and 1/2, respectively. 

The inclusive electron scattering cross section is related via the optical theorem 
to the imaginary part of the forward doubly virtual Compton amplitude 7*A^ — > 
7*A'^. Starting from a generalization of Equation E31 for virtual photons, Drechsel 
et al. defined dispersion relations between the (e, e') partial cross sections and 
forward doubly virtual polarizabilities p45j) 

o-t(j^,Q )^ 



[a + l3]{Q\Q') = 



lA{Q^)aQED/M^ 

lo{Q\Q^) 



1^0 



1 

AMuqed 

1 

2^ 

16M3 



2xFi{x,Q'^)dx 



OLQED 



dx 



|4M2 



[F2 - 2xFi] + x^F2 



1 



(2vr)2 
1 

2^ 



K du 

— (JTT — 



(73) 

(74) 
(75) 



2^0 



IGM^aQED 



1 1 

27r2 



dxx 



X 



gi{x,Q ) g2{x,Q ) 



o-Lr(i^, Q' 



1^0 



T 

dv 
> — 



(76) 
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= ^-^^^^^Jdxx'[g,{x,Q')+g,{x,Q')] (77) 



A discussion of the Generalized GDH sum rule I a (Equation [75]) is beyond the 
scope of this review. Experimental determination of these forward polarizabilities 
requires a full separation of the partial cross sections in unpolarized and polarized 
lepton scattering. However, the polarizability a + f3 dominates the unpolarized 
cross section, and at high Q^, aL/[a + P] oc illustrating that a^, is a 

twist-four matrix element (|145|) . The ^ limits of a + /3 and 70 are the 
usual Baldin sum rule (Eouation I26() and the forward spin polarizability relation 
(Equation 155)) . Using the Wandzura & Wilczek (|217|) estimate of §2 ■ 

hT{Q\Q^) ^ 7o((3',(3')/3, asQ^^oo, (78) 

the Jefferson Lab Hall C collaboration has separated the Fi and F2 structure 
functions on the proton by the standard Rosenbluth procedure of measurements 
at fixed and variable e (j224i|l . Figure El displays the results for the forward 
polarizability [a + /3](Q^Q^). The MAID calculation (I145U157I) . with the vriV, 
rjN , and vrvrA^ intermediate states included is in good agreement with the data, 
if the integration is truncated to W < 2 GeV. The estimate from DIS structure 
functions is in good agreement with the data if the integrand in both cases is 
truncated to > 2 GeV. The Jefferson Lab Hall B and C collaborations have 
also measured the moments of the F2{x,Q'^) structure function, using previous 
estimates of the ratio R = ai/aT (|2181I219|1 . 

The Jefferson Lab Hall A GDH collaboration has separated the and g2 struc- 
ture functions from measurements of the ^He(e, e') reaction (|22H1 . Figure ITSl 
displays the forward spin polarizabilities of the neutron. The MAID curves 
under-predict the data for 70 at low Q^, but are otherwise in good agreement 
with the data. The curves from Kao et al. are 0(j>^) HBChPT ()222j) . The curves 
of Bernard et al. are 0{p^) Relativistic Baryon ChPT, with the shaded bands 
including the effects of explicit inclusion of the A and vector mesons ()223j) . These 
explicit resonance effects are large for 7q , and bring the ChPT predictions into 
agreement with the data at = 0.1 GeV^ (but already diverge from the data 
by = 0.25 GeV^). Although the effect of the resonances is much smaller for 
Slt, the ChPT predictions disagree sharply with the data in Figure ITHl However, 
the disagreement is comparable to the differences between the two 0{p'^) ChPT 
calculations, and the change from Oi^p"^) to 0{p^) ChPT calculations H222|) . 



4 SUMMARY, OUTLOOK AND CONCLUSIONS 

Recent advances in polarized electron sources, polarized nucleon targets and 
nucleon recoil polarimeters have enabled accurate measurements of the spin- 
dependent elastic electron-nucleon cross section. New data on nucleon electro- 
magnetic form factors with unprecedented precision have (and will continue to) 
become available in an ever increasing Q^-domain. The two magnetic form fac- 
tors and G^j closely follow the simple dipole form factor Gd- G^e/^\i drops 
linearly with and appears to drop at the same rate as from ~ 1 GeV^ 
onwards. The Q^-behavior of has provided a signal of substantial non-zero 
orbital angular momentum in the proton. Only scant data are available in the 
time-like region. The full EMFF data set forms tight constraints on models of 
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Figure 17: Q^-dependence of the generalized forward Baldin sum rule (Equa- 
tion EHI- The real photon point is from (|14H1 . the Jefferson Lab data are from 
(|224)1 . The curves are taken from Figure 6 of (|145)1 . The sohd curve is from the 
MAID parameterization of 7(*)A^ ttN multipoles (|157jl . The dashed curve in- 
cludes the contributions of ijN and vrvrA^ intermediate states (|145|) . In the lower 
panel, the dotted curve is the estimate from the DIS structure function Fi (,220,1 , 
for > 2 GeV. 



nucleon structure. So far, all available theories are at least to some extent effec- 
tive (or parametrizations). Still, only few of these adequately describe all four 
EMFFs. Only lattice gauge theory can provide a truely ab initio calculation, but 
accurate lattice QCD results for the EMFFs are still several years away. A scaling 
prediction has been developed for the ratio of the Pauli and Dirac form factors, 
which the data appear to follow even at a Q^-value as low as 1 GeV^. Novel 
procedures allow a visualization of the nucleon structure as a function of the 
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Figure 18: Q^-dependence of the neutron generalized forward spin polariz- 
ability 7o(Q^)Q^) (Equation [7^ and fongitudinal-transverse spin polarizability 
S^lxiQ"^ 7 Q"^) (Equation [771) . The shaded bands on the the systematic 

error bands of the Jefferson Lab data ()221. 225). The (dark) sohd curves ex- 
tending beyond =1 GeV^ are the MAID parameterization of the 7*A^ — > vrA^ 
ampUtudes H145LI157)) . The other curves (including the shaded bands) are ChPT 
calculations (see text for details) H2221I223|) . 

momentum of the struck quark. A fully three-dimensional picture of the nucleon 
will become available when future exclusive data have allowed the determination 
of the Generalized Parton Distributions. 

Measurements of the nucleon polarizabilities have followed a 40-year odyssey that 
parallels the history of form- factor measurements. Powerful new theoretical and 
experimental techniques allow high precision measurements of the polarizablities 
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and their Q^-dependent generalizations. These results offer new evidence for the 
interplay of constituent quark and pion degrees of freedom at modest distance 
scales within the nucleon, and the dynamics of the elementary current quarks 
of QCD at large Q^. This program will continue with both existing and future 
experimental facilities. 
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